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We present a model example of a quantum critical behavior of renormalized single-particle Wan- 
nier function composed of Slater s-orbitals and represented in an adjustable Gaussian STO-7G basis, 
which is calculated for cubic lattices in the Gutzwiller correlated state near the metal-insulator tran- 



o 

I sition (MIT). The discussion is carried out within the extended Hubbard model and the method 

, of approach proposed earlier [cf. Eur. Phys. J. B 66, 385 (2008)]. The component atomic-wave- 

^ ' function size, the Wannier function maximum, as well as the system energy, all scale with the 

increasing lattice parameter R as [{R — Rc)/RcY with s in the interval [0.9, 1.0]. Such scaling law is 
interpreted as evidence of a dominant role of the interparticle Coulomb repulsion, which for R> Rc 
■ is of intersite character. Relation of the insulator-metal transition lattice-parameter value R = Rc 

to the original Mott criterion is also obtained. The method feasibility is tested by comparing our 
results with the exact approach for the Hubbard chain, for which the Mott-Hubbard transition is 
absent. In view of unique features of our results, an extensive discussion in qualitative terms is also 
provided. 



C/2 , PACS numbers: 71.30.-l-h,71.28.-l-d,71.10.Fd 



1. Introduction 



i Metal-insulator transition (MIT) of the Mott-Hubbard type represents one of the central problems in condensed 
Q \ matter physics and materials science, since it exemplifies a transition from well defined atomic (localized or confined) 
O . states to delocalized (Bloch- or Fermi-liquid) type of states in a solid^ and in other systems.— The simplest model 
system represents a lattice of spin-1/2 fcrmions with one particle per atomic site. Mott expressed this localized- 
^SJ . delocalized transformation in terms of critical density nc of fcrmions (or equivalently, in terms of critical interparticle 
\ distance Rc — n^^^"^ , for three-dimensional system and the effective Bohr- radius a of the atomic states formed at the 

' transition, in the form ri^^^^a ~ 0.2 — 0.25. On the other hand, Hubbard- and others'* formulated the criterion in 
, the form that at the transition the magnitude of the intraatomic Coulomb repulsive energy among the particles (L/) 
' is equal to their kinetic (band) energy characterized by the bare bandwidth W , i.e. that U ~ W . These two criteria 
i are regarded as related (cf. Appendix A). 

The metal-insulator transition is also well defined as a phase transformation in thermodynamic sense with the 
^\ \ corresponding critical points on high- and low-temperature sides»^ As close to MIT the band energy (negative) 
i is almost compensated by the repulsive (Coulomb) interaction, the much smaller thermal and for atomic-disorder 
energies can drive the metallic system towards the state with localized spins. In effect, the resultant phase diagrams 
contain first-order pressure-temperature transition line with a classical critical point at temperature Tcr > and a 
quantum critical point at T = 0. The existence of the former critical points was confirmed experimentally recently^ 
and shown to represent properties of the Van der Waals liquid-gas critical point for the three-dimensional case of 
V203:Cr. The basic question remains under what conditions a quantum critical point (QCP) appears also at Tcr = 
between the metallic and insulating phases, as its existence marks explicitly the fundamental boundary between the 
atomic and the condensed (delocalized, deconfined) quantum states. The appearance of such QCP is usually obscured 
by the presence, at low temperatures, of antiferromagnetismi^ Nonetheless, the answer to the above question, even 
without inclusion of magnetic ordering, would delineate the basic characteristics of the (quantum) critical point at 
the border of localization driven by the interparticle interaction.^ 

MIT involves the drastic change of quantum macro states from a metal to an insulator, so a natural and not 
addressed, as far as we are aware of, is the question how this drastic change of the macro state is refiected in the cor- 
responding single-particle wave-function change when approaching the transition from either side. The parameterized- 
model approaches^i^ leave the Wannier functions determining the microscopic parameters such as the hopping integral 
t or Hubbard interaction U, as fixed when crossing the transition. On the other hand, the original Mott approach^ 
introduces intuitively the concept of an emergent atomic state at the instability of the metallic state (cf. Appendix 
A). Thus, the missing question is: how the electron correlations and the single-particle states are interrelated micro- 
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scopicallyl This question is particularly important, since the Fermi-liquid state near the Mott-Hubbard localization 
is particularly robust and therefore the Mott-Hubbard localization At T > has a strongly discontinuous nature^, 
apart from the classical critical point of Tc- the crossover behavior for the temperatures above it. Furthermore, at 
the transition the renormalized band and the Coulomb energies are comparable and compensate each other?^ 

In this paper we provide a simple answer to the above question by treating electronic states close to MIT within 
the approach in which both the interelectronic correlations and the single-particle Wannier wave functions {wi{r)} 
are treated on the same footing^. By this we mean that those wave functions are determined simultaneously with 
the interparticle correlations determining the many-body ground state. As a result, we obtain a singular behavior 
of the single-particle (Wannier) wave function renormalized by the correlations. This nonanalytic character, and 
concomitant with it unique scaling laws, appear only when the wave function is determined explicitly a posteriori, 
i.e. is readjusted in the correlated state. Such a simultaneous determination of the single-particle wave function 
and two-particle aspects of the many-particle dynamics is indispensable in the situation when the single-particle and 
interaction is comparable or, in some cases, the latter even becomes dominant. In this sense, our analysis complements 
that performed within either LDA+DMFT^° or LDA+\}^ methods, as well as the QCP studj*^, all of which contain 
parameters characterizing electronic correlations introduced only after the LDA or similar calculations have been 
carried out. Obviously, the Gutzwiller- ansate used here provides an exact wave function only in the limit of high 
dimensions (i.e. above the upper critical dimension, not determined as yet), but it provides at least a mean-field- 
like discussion of the renormalization aspects of the single-particle wave function near MIT, not discussed so far in 
the literature. In connection with this last limitation of our method, we would like to point out, that it would be 
very interesting to apply methods such as DMFT^*' to the proposed optimization of the Wannier functions in the 
correlated state. However, no reliable method of evaluating explicitly the DMFT ground-state energy as a function 
of microscopic parameters has been proposed so far, starting from which the explicit wave-function optimization 
can be undertaken. We should also underline, that our method incorporating an exact diagonalization was applied 
earlier to both correlated nano^ and macro— systems. We should also point out that this method is free from the 
double counting of the repulsive Coulomb interaction. Therefore, even though the presented below results describe a 
model situation and within a rather simple scheme, they can be regarded as an essential additional ingredient to be 
implemented in modeling real correlated systems near MIT. As was underlined earlier^^, within the present method 
we can follow the system evolution as a function of interatomic distance, not only as a function of model parameters. 
Such circumstance allows for an explicit analysis, at least in a qualitative manner, of the system evolution as a function 
of external pressure. 

The structure of this paper is as follows. In Section 2 we introduce briefly the method devised before i^ i^^'^'^ In Sec. 
3 we introduce an original critical scaling of the single-particle wave function characteristics and of the ground-state 
energy. In Sec. 4 we show that the critical scaling is absent for one- and two-dimensional lattices, in the former case in 
agreement with the absence of metal-insulator transitions, as discussed rigorously by Lieb and WuM^ In this Section 
we also derive the classic Mott criterion of the metal-insulator transition within the frame of our approach. Sec. 5 
contains both a brief discussion of extension to more realistic situation and concluding remarks. In Appendix A we 
provide an elementary derivation and discussion of the original Mott criterion.-, whereas in Appendix B we define an 
adjustable Gaussian basis in which the renormalized Wannier function has been determined explicitly. 

A methodological remark is in place here. In our earlier work,^^'^^ the method of approach composed of the 
diagonalization in the Fock space combined with the single-particle wave-function renormalization in the Hilbert 
space, has been applied to the discussion of the metallic state up to the metal-insulator transition point. No scaling 
presented here was either noticed or discussed thereJ^ Therefore, we demonstrate here, in an explicit manner, the 
applicability of our method when crossing the metal-insulator boundary of the Mott-Hubbard type. Namely, we 
study explicitly the gradual evolution of renormalized Wannier functions into their atomic correspondant in the large 
interatomic-distance limit. On the two above aspects rests our motivation that the present paper extends essentially 
our earlier approach i^J^ii^ 



We start from extended Hubbard model with inclusion of intersite Coulomb interaction and ion-ion interaction as 
represented by the parameterized Hamiltonian^-^^'^'^ 



2. The method: adjustment of the wave function in the correlated state 




(1) 
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where the first term describes the effective atomic energy with 

= + ^ ^ (i^., + 2/i?,,) , (2) 

the second represents the hopping between the nearest-neighboring sites, the third the Hubbard (intraatomic) and 
the fourth is a part of intersite Coulomb interaction, with Srii = 1 — rn being the deviation from the integer electron 
occupancy = 1, and Kij - the intersite Coulomb interaction. The term {2/Rij) in ^ expresses the classical 
Coulomb repulsion (in atomic units) for two ions separated by the distance Rij. It is important to note that the 
presence of intersite Coulomb interaction is necessary in achieving a proper atomic limit when i? — > oo. The model 
has been defined in details earlier.— i^^'^'^ 

In the Mott-Hubbard insulating ground-state the hole-hole correlations are absent i.e. ((1 — ni){l — rij)) = 0. The 
role of these correlations is indeed negligible if the fundamental correlation function d = {nii-^nii) vanishes or is very 
small, as also discussed earlier It vanishes in the Mott insulating state (for R > Rc). 

The microscopic parameters of this model are expressed via the Wannier functions in the standard manner— i^^, 
{wi{r)} = {w{r — Ri)}, as follows: the bare atomic energy is = {wi\Hi\wi) , the hopping integral tij = (wilHilwj) , 
intraatomic-interaction magnitude U = {wiWi\Vi2\wiWi) , and interatomic-interaction magnitude Kij = {wiWj\Vi2\ 
WiWj), where Hi is the Hamiltonian for a single particle in the system, and V12 represents interparticle interaction. The 
Wannier functions are expressed in terms of adjustable Slater atomic functions, i.e. Wi{r) = /3\I'i(r) — 'yJ2j=i 

where z is the number of nearest neighbors, P and 7 are mixing coefficients, and = (a^/vr)^/^ 

ea;p(— a|r — Ril) is the Is Slater function centered on the site i = R^. In the concrete calculations, they are represented 
by adjustable Gaussians (e.g. ST0-7G basis in dimensions I? = 1, 2, and 3, cf. Appendix B). The properties of those 
functions have been discussed in detail beforei^ii^ Here we concentrate on the novel aspects of the approach only. This 
means we explicitly determine the critical interatomic distance R = Rc for the metal-insulator transition, as well as 
determine the wave function characteristic as a function of R, both above and below Rc- Based on these results we 
determine scaling properties of physical quantities vs. the parameter (i? — Rc)/Rc- In effect, we show that not only 
the wave-functions determination and the electronic correlations are interrelated, but also the former characteristic 
exhibits a singular behavior at the onset of Mott-Hubbard localization. This last feature of our results represents the 
principal message of this paper. 

3. Quantum critical behavior of the wave function 

3.1 Overall features of the wave function: evolution of Wannier functions into atomic wave function with 

i? — > 00 

As said above, the fundamental principle behind our approach is that the wave functions {wi{r)} determined 
variationally in the Hilbert space are treated on the same footing as the diagonalization of ^ in the Fock space. Such 
diagonalization is possible in an exact manner for nanoscopic and infinite Hubbard chain {D = 1) systems onlyJ^ii^ 
However, also the Gutzwiller wave-function (GWF) and the Gutzwiller-ansatz (GA) approximations lead to close 
results in the latter case (see below) provided the wave functions {^^(r)} are properly readjusted in the correlated 
state to achieve the ground-state energy as a global minimum for given lattice parameter R. In other words, the energy 
of the correlated state is readjusted iteratively multiple times by adjusting the component atomic-wave-function size 
and in effect, also the microscopic parameters £eff, tij, U, and Kij determining the many-particle ground state. 
As a result of such iterative approach, we obtain renormalized wave functions {^^(r)} and the microscopic parameters, 
as well as the ground-state energy, and most importantly, an explicit scaling of the physical properties, all as a function 
of R. Below we analyze first the results obtained for GA for three-dimensional cubic lattices, before testing the method 
feasibility for _D = 1 and 2 situations. The main emphasis is lead on the novel singular scaling properties near the 
Mott-Hubbard critical spacing R = Rc. 

In Fig. [T] we plot exemplary Wannier function centered at R^ = for simple cubic (SC) lattice and drawn along 
[1,0,0] direction and for R < Rc (dot-dashed line) R — Rc (solid line), and R ^ Rc (dashed line). One observes 
immediately a surprising feature, namely the wave function for R > Rc is more extended than that for R — Rc and 
with the increasing R the size a — 1/a gradually reaches its Bohr atomic Is-state value oq. An overall wave- function 
profile evolution as a function of i? is displayed in Fig. [51 On the whole, the wave function seems to evolve continuously 
for R > Rc and reaches its atomic Is shape gradually, except for a local maximum which appears at the specific lattice 
constant R = Rc, as marked by the vertical arrow. Before studying this unusual feature let us note, that our method 
enables us to study this wave-function evolution in the whole range of R (across the metal-insulator boundary) as 
long as the tight-binding approximation is reliable. This is so, as we shall see and the interesting scaling laws appear 
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FIG. 1: (Color online) Exemplary shapes of the Wannier function wo(r)centered at the site Ri = for simple cubic (SC) 
lattice along [1,0,0] direction and for three lattice spacings marked: R < Rc (dot-dashed line), R — Rc (solid line), and R > Rc 
(dashed line), ao is the Is Bohr radius. 

near Rc which falls in the limits of large and small interatomic separation, R 3> and R < a^^, respectively. 
The value R = Rc is determined from the point when the renormalized band energy is completely compensated by 
the Coulomb interaction part and coincides with the correlation function d =< ni^riii in the Gutzwiller-ansatz 
approximation for the half-filled-band case. 

3.2 Wave-function singular behavior at Mott point and scaling properties in the critical regime R Rc 

To study the noted above nontrivial behavior of the wave function in detail, we plot in Fig. [31 the relative inverse 
size 5a/ a = \a{R) — a{Rc)\/a{Rc) of the component atomic functions as a function of relative lattice spacing 5R/ R = 
[R — Rc)/Rc in the regime R ^ Rc- A clear da/dR discontinuity (cf. inset) is observed at R = Rc for all cubic 
lattices, a rather unique and unexpected feature, which is completely absent for the case with bare (unrenormalized) 
wave functions, as then a~~^{R) would be simply a constant independent of R. Note also that the wave function is the 
narrowest at the Mott-Hubbard transition, at which d =< nii^riii (d > in the metallic phase). To determine 
the universality of the behavior, we have replotted in Fig. 2^ the result of Fig. [3] for R > Rc in a doubly logarithmic 
scale. One observes a clear power law scaling 5a/ac ~ [{R — Rc)/RcY, with s — 0.96±0.01 for not too large R above 
Rc- The results for R < Rc exhibit also a similar trend with the exponent s with an overall exponent s slightly lower 
(s ~ 0.93) for SC and BCC lattices i.e. weakly dependent on the type of cubic lattice selected, as shown in Fig. [3Jd. 
This nonuniversal character (cf. also the part for R < Rein Fig. [3]) may be caused by the circumstance that the tight- 
binding approximation for the wave function should be extended to more distant neighbors, as the presence of carriers 
screens out effectively the Coulomb interaction (the value d = {ni-^riii) increases as R decreases). Nonetheless, in spite 
the lack of simplicity of the precise scaling for i? < i?c, it is reassuring that the critical-exponent values for R < Rc 
are close to that for R > Rc- This power-law scaling describes the singular behavior of the aioTnic-wave- function size 
a = 1/a. Also, the function in Fig. [3] with R oo does not reach zero, since a{R — > oo) = Gq^ > a{R = Rc). 

One can also study the related properties of the renormalized Wannier function directly. For that purpose we have 
displayed in Fig. [5] the maximal value ofwi{r = 0) as a function of the scaling parameter {R — Rc)/ Rc- Again, also for 
that quantity, we observe a similar type of scaling as for a~^, with the critical exponent a = 0.92 ± 0.01 (cf. inset). 
However, the corresponding scaling law in the regime R < Rc is again not as clear, as that for R > Rc (cf. Figs. [3] 
and ini). Note that close to R ^ Rc we have (a — ac)/<^c — (a — ac)/ac, where a is the renormalized Bohr radius in 
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FIG. 2: Overall space profiles of the renormalized Wannier function for simple-cubic lattice, as a function of lattice parameter 
R. The arrows mark the position of the singular (critical) behavior elaborated below in main text (cf. also Fig. |3l The Wannier 
function evolves with i? — > cx) (i? S> Rc) into atomic Is wave function. 
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FIG. 3: (Color online) Relative inverse wave function size 5a/ac as a function of relative lattice parameter 5R/Rc for simple 
cubic (SC), body centered cubic (BCC), and face centered cubic(FCC) lattices. The Mott-Hubbard transition is marked by the 
vertical dashed line. Inset: derivative da/dR vs. 5R/R with a singularity at SR = R — Rc = 0. Note the universal behavior in 
the range R> Rc- 
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FIG. 4: (Color online) a) Detailed scaling 5a/a vs. 5R/R and in the doubly logarithmic scale for R > Rc (top) and for R < Rc 
(bottom). The straight lines represent the fitted curves ~ (R/Rc)". b) The scaling for R < Rc- Note absence of the degree of 
universality in the alter (b) case. 

the correlated state. In other words both Sa/ac and Sac scale the same way with SR/ Re- 
There is no obvious quantity like a^^ characterizing the width of the renormalized Wannier function Wi{r). This 
is the reason why we calculated an overaU (average) size of the renormalized Wannier function defined as (r) = 
J rw'^{r)d^r, and have plotted it in Fig. [S] as a function of SR/Rc for the three cubic lattices. While we observe 
anticusp behavior of (r) at R = Rc for SC and BCC lattices, there is no clear sign of that happening for the FCC 
case. However, one has to realize, that calculation of (r) involves a three-dimensional integration over space in the 
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FIG. 5: (Color online) Maximum iuq (0,0,0) of the Wannier function vs. SR/Rc- Inset: detailed scaling of the renormalized 
Wannier function maximum in the doubly logarithmic scale (of same type as in Fig. |4]). Fitted straight line (for R > Rc) 
represents the function ~ [{R — Rc) / Rc]'^'^^ ■ Note the essential similarity of the scaling above and that shown in Fig[3]and|31 

situation when the Wannier functions are strongly anisotropic— Because of this specific reason, we regard the scaling 
of wo(r) maximum shown in Fig. [5] as more direct, as since involves the wave function characteristic which is not 
averaged out and, obviously, independent of the spatial direction. Nonetheless, in Fig. [7]wc display the corresponding 
scaling of the relative difference ((r) — {vc))/rc for R > Rc- The less systematic behavior even for R > Rc does not 
allow us to draw any definite conclusions for that regime. Note that for R > Rc even in the case of FCC lattice, where 
(r) is approximately continuous when crossing R^ we observe the same type of scaling, as that discussed for SC and 
BCC lattices. This means again that the scaling in the Mott insulating phase is induced by the intersite Coulomb 
interaction. 



The scaling with varying 5R/R is also observed for the ground state energy, as shown in Fig. [5] The approximate 
dependence ~ {SR/RY of Eq when approaching the atomic limit can be attributed to the dominant role of the 
Coulomb repulsion (note the absolute value of the energy). 

To summarize this Section, both the Figs [4] and [7] demonstrate (for R > Rc), respectively the appearance of 
a overall scaling of the i?^" type, with n ~ 1. This dependence could be seen in the pressure dependence of the 
orbital size when close to the critical value R — Rc- However, the relative changes of a and (r) are rather subtle and 
the question remains if they can become observable in the present day e.g. scanning tunneling observations of the 
localized orbitals in the Mott insulating state. For such purpose, relevant are the plots of the maximal value wo{0) 
of the Wannier function relative to that aX R — Rc (cf. Fig. [5]). The particle occupancy change at r = reaches up 
to 20% of the peak value upon change of 20 — 25% of the lattice constant either way, so it could become observable. 
On the basis of our results one can also say, that our approach provides evidence for effects in the Wannier function, 
which are pronounced already in the Mott-Hubbard insulating state. In other words, the Hubbard split subband 
picture of the Mott insulator is concomitant with a strong renormalization of the wave function characteristics when 
approaching the transition point. This critical behavior does not alter the monotonic increase of the U/\t\ ratio with 
the increasing lattice constant*i2,, as shown explicitly in the Fig. [9] for all three cubic lattices. The almost exponential 
growth of this ratio for R > iag is a signature of entering into the strong-correlation limit. Also, the presence of the 
scaling demonstrates the onset of the collective character of the quantum critical regime, in which the Mott insulating 



3.3 Scaling property of ground-state energy and physical meaning of the scaling 
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FIG. 6: (Color online) First moment (an average size) of the renormalized Wannier function for cubic lattices. Note the same 
type of discontinuity as in Fig. |3]for SC and BCC lattices. 
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FIG. 7; (Color online) Relative average size ((r) — {f'c))/ {tc) of the renormalized Wannier function for SC, BCC, and FCC 
cubic lattices as a function of relative lattice constant SR/Rc. The discontinuities appear at the critical value Rc. The — > oo 
asymptotic value of (r) = 1.5ao is that for atomic Is wave function, with ao being the Bohr radius. 
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FIG. 8: (Color online) Same scaling as in Fig. |4] but for the ground state energy; the fitted exponents are listed for the three 
cubic lattices. In the asymptotic limit SR/Rc >> 1 the three values of relative energy merge, since the system approaches 
the atomic limit. Close to the atomic limit the scaling is not obeyed. This means that a collective Mott state sets in with 
R — > -Rc^ . 



state cannot be regarded as a collection of atomic states with unpaired spins, at least when R ^ R^. 



4. Effect of lattice dimensionality: feasibility of the method and the Gutzwiller-ansatz approximation vs. the 

original Mott criterion 

Finally, we examine the role of lattice dimensionality on the singular behavior of the wave function. For that purpose, 
we plot in Fig. [To] the inverse wave-function size a ^ vs. R/ao for the linear chain within the modified exact Lieb-Wu 
(LW)fi^ Gutzwiller- wave-function (GWF),^'^ and Gutzwiller-ansatz (GA)i^ii^ii^ solutions, respectively. None of the 
three methods provides in D = 1 case (and GA also in the D — 2 case) any type of singular behavior observed in 
Figs. [3] and [5] for the D = 3 case. What is more important, all the three methods provide quite similar results for 
D — 1, i.e. there is no critical behavior of the type discussed abovei^ This circumstance provides us with some 
confidence, that modified by us GA method (including the interatomic interaction exactly in the Mott insulating 
state) bears some physical relevance to the MIT transition treated in mean-field approximation. MIT transition for 
D — 1 and D — 2 lattices which appears in the standard (mean-field) Gutzwiller ansatz, does not show up in the 
behavior of the renormalized wave function in real space. This is one of the main supplementary results of this paper. 
To illustrate that, we plot in Fig. [TO] exemplary Wannier function vs. R/uq for triangular (TR) and square (SQ) 
lattices. The results for D — 2 are obtained for GA solution. Again, no critical behavior discussed above is observed. 

On the basis of these results we draw the conclusion, that within our approach the critical behavior of the Wannier 
function appears only for high-dimension (D > 2) lattices. This circumstance is in accord with the notion that GA 
results are more realistic in high dimensions. The striking coincidence of these results with the critical dimensionality 
{D > 2) for the onset of Landau- Fermi liquid stability should therefore be mentioned.-^ The metallic 3D state of the 
present approach is represented by an almost localized Fermi liquid^ (ALFL), which differs from the Landau Fermi 
liquid by a strong (and spin dependent) mass enhancement factor, as well as by a presence of metamagnetism.^® 

As a supplementary information, we provide in Fig. [TO] the average size (r) of the renormalized Wannier function 
for linear chain (CH) and the two-dimensional lattices, SQ and TR. As is the case in Fig.[Tll no singular character has 
been observed for this quantity. These results confirm again that our results are sensitive to the lattice dimensionality 
and nontrivial properties appear only in the high-dimension limit, where the mean-field (Gutzwiller) approximation 
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FIG. 9: (Color online) The U/\t\ ratio as a function of the interatomic spacing for cubic lattices showing its monotonic increase 
upon increasing R, unlike the scaling of the Wannier function characteristics. 




FIG. 10: (Color online) Inverse size of the renormalized atomic wave function for linear chain (CH) vs. R/ao within the 
exact Lieb-Wu (LW), the Gutzwiller- wave-function (GWF), and the Gutzwiller-ansatz (GA) solutions. No singular behavior 
is observed at any R. Inset: Wannier-function maximum it)o(O) vs. R/ao for linear chain (CH), square (SQ), and triangular 
lattices (TR), with no discontinuity detected in either of the cases. The scaling in GA approximation discussed above does not 
apply here. 
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FIG. 11: (Color online) Space profiles of the Wannier functions for triangular (TR, top) and square (SQ, bottom) planar lattices, 
both as a function of lattice constant R. No singular behavior within GA approach is observed. Note a gradual evolution with 
the increasing of R of the Wannier functions into the corresponding atomic (Is) form. 

may be regarded as at least qualitatively correct.— Here they appear already for D = 'd. 

One additional basic feature of our approach should be mentioned. Namely, in Table I we list characteristics of 
the Mott-Hubbard transition for the cubic lattices. These results intercorrelate nicely with the Mott criterion in the 
following manner. The critical carrier concentration for the Mott transition in the present situations is nc — ^/Rc 
for SC, 2/Rl for BCC, and 4/i?g for FCC. Thus nc = n/R^, with n = 1,2, and 4, respectively. The Mott criterion 

discussed in qualitative terms in Appendix A, takes for those lattices the form: n^"^ ■ a = v}^^ / (Rcac) = 0.20, 0.26, 
and 0.33, for SC, BCC, and FCC lattices, very close values to that obtained originally by Mott,^ and rederived in 
an elementary manner in Appendix A. This connection provides an additional argument for a qualitative correctness 
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FIG. 12: (Color online) Average size of the renormalized Wannier function for one- (CH), and two-dimensional (SQ, TR) 
lattices. No singular behavior observed apart from a slight slope change for the linear chain. 



TABLE I: Microscopic parameters near the mean-field critical point calculated for the cubic lattices. 



Struct. 


aao 


Rc/ao 


(U/W)c 


SC 


1.099 


4.236 


1.337 


BCC 


1.109 


4.384 


1.080 


FCC 


1.128 


4.351 


0.880 



of our approach in the sense, that a weU defined value of obtained from GA solution, corresponds nicely with 
the historic criterion introduces by Mott.- Parenthetically, the agreement (cf. Table 1) relates nicely the Mott and 
Hubbard criteria for MIT. 



5 Outlook 

5.1 A brief overview of the method combing first- and second-quantization schemes 

Our whole method starts with the definition of the field operators 

r $(r)=i:*»(r)(a.T,«.i), 
|*(r)=E<f.(r)(^ 

and of the corresponding Hamiltonian (1) in the second quantization representation. The method of approach is 
equivalent to that starting from N-particle Schrodinger equation provided the single-particle basis {$i(r)} is complete 
in the quantum-mechanical sense. The basis {$i(r)} can be arbitrary, even not orthonormal^^, although the choice of 
orthogonal basis simplifies the whole procedure. Within this scheme, the starting Hamiltonian (1) represents a model, 
as it includes only a single Is Wannier function {^^(r)} per atomic site (no other orbitals included).. 
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The question is how to select the correct wave- function basis Wi (r) for this model. This question has been addressed 
a number of times in a general case at the original stage of setting the quantum theory.^^ In general, by introducing 
the occupation-number representation we determine combinatorially (or algebraically) the number of single-particle 
configurations and the resulting interaction energy as an average over the set of these coherent configurations (a 
multiconfigurational average). The single-particle wave function on other hand, should be adopted variationally in the 
resultant multiparticle state, particularly if the interaction energy is comparable to the single-particle (band) part, as 
no perturbational approach starting solely from the single-particle part is appropriate. Conversely, the single-particle 
wave-function size readjusted in the correlated state represents an additional factor lowering the system energy. ■^■^ This 
procedure has been carried out successfully first in the case of the He atom and H2 moleculci^ Minimally, determining 
the system evolution as a function of model parameters, which is insufficient from the physical point of view. Obviously, 
such a variational wave-function readjustment is carried out already in e.g. LDA method^^ Our method of approach 
demonstrates explicitly, albeit in a model situation, that the single-particle wave function variational readjustment 
only is insufficient, as an intrinsically different ingredient, the two-particle correlation function, is ignored in the 
process. The LDA and related methods are certainly valid (and often sufficient) when the two-particle correlations 
can be approximated by a product containing density of particles only, i.e. in the weakly or moderately regime of 
correlated particles. 



In summary, we have extended the standard treatment of the Mott-Hubbard transition by incorporating into it a self- 
consistent scheme the single-particle wave-function renormalization, as well as demonstrating its singular behavior for 
three-dimensional lattices within the Gutzwiller approximation. These effects are absent in lower dimensions (D < 3). 
We have also related the critical interatomic distance Rc to the original Mott criterion. The correlation effects in 
the single-particle wave-function-shape develop already in the Mott-Hubbard insulating phase. The wave-function 
characteristic length a^^ (or its average size (r)) play the role of a coherence length for a single- fermion placed in 
the mean field of all other particles. In connection with this one has to mention the works on the scaling theory of 
Anderson transition, where the coherence length scales with the magnitude of interaction with the exponent v = —1. 
In our case roughly ((r) — (rc)) ~ (i? — Rc) ~ {Kc — K) / for R > Rc, where Kc is the value of K for R = Rc- This 
means the scaling is different in the present case and differentiates between the Mott and Anderson transitions. 

The present approach can also be extended to the orbitally degenerate 3d states along the lines proposed earlier J^i^ 
Important is the circumstance that the present method, although still at a simple modeling stage, avoids the notoriety 
of counting twice the interparticle interaction, present in any of the current ab initio approach, and goes beyond the 
parameterized-model calculation o^^i^^ , as it provides the system evolution as a function of lattice constant. Also, it 
is of fundamental importance to try to extend the present hypothesis of a quantum critical scaling to check if it holds 
beyond a mean- field (high-dimensional) treatment of strongly correlated fermions near the Mott-Hubbard transitions. 



The authors are grateful to Krzysztof Byczuk for interesting points on DMFT method. The work was supported 
by the Grant No. NN 202 128 736 from Ministry of Science and Higher Education. Also, the research is performed 
under the auspices of the ESF Network INTELBIOMAT and the National Network Materials with Strongly Correlated 
Electrons. The authors acknowledge also discussions with Dr. Adam Rycerz on and adjustable Gaussian basis. We 
are also grateful to Dr. Dominique Delande from Universite Paris VI for turning our attention to Ref.— . 



The Mott criterion of electron localization can be visualized in an elementary manner from the following reasoning 
of the gas instability. The kinetic energy per particle in 31? ideal gas is 



5.2 Conclusions 



Acknowledgments 



Appendix A: Elementary meaning of the original Mott criterion 




(3) 
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where is the gas Fermi energy p = N/V-the particle density and m*-the effective mass in the gas. The electron- 
electron interaction at the onset of localization can be estimated as 

where the classical interparticle distance has been taken as de_e — {N/Vy/"^ — p^^l'^ . 

Now, the gas instability point is defined for a critical density p — pc for which e = (e-e- This condition leads to 
identity 

-r^e =o ^^0.17, 5 



In other word 



aepl'^ 0.2 (6) 



which represents the Mott-Wigner criterion of localization, with as being the effective Bohr radius of Is-type of 
atomic state for electrons in the medium with dielectric constant e and with effective mass m*. For p < pc the 
Coulomb repulsion dominates (e < Ce-e)- 



Appendix B: Approximation of Slater Is functions by an adjustable STO-nG basis 

In our previous paper— regarded as Part I we have not defined explicitly the adjustable Gaussian basis. The details 
of this are provided below. 

The method of approximating the molecular orbitals by Gaussians was introduced by Hehre at al.-^. In this 
representation the Gaussian orbitals have been introduced in the form 

5l'^^(r)=(^)'''e-"l^---^ (7) 

where is the reference atomic site and the index "a" labels different functions. In the basis STO-nG the Is wave 
function is written as a linear combination (contraction) of n such Gaussians. As we would like to have a variable-size 
orbitals, we make a readjustment Fa — » rjj = aVa- In effect, the atomic Is wave function has the form 

n 
a=l 

The parameters Pa and T'^ are determined from minimalization of energy of a single atom described by the following 
Hamiltonian in the atomic units: 

and are independent of a. In Ref;^^ one can find the expressions for various integrals involving Gaussian basis. 
Particularly important for us are the following integrals: 



1. The overlap integral between sites i and j 



2. The kinetic energy matrix elements 
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TABLE II: Values of parameters /3a and Fa for a = 1, 7 for Is Slater wave functfon obtained in (STO-nG) representation for 
n = 3, 5, and 7, respectively. 





STO-3G 


ST0-5G 


STO-7G 


a 


Pa 


i a 


Pa 


i a 


Pa 


i a 


1 


0.7079 


0.4037 


0.4862 


0.3429 


0.3348 


0.3073 


2 


0.3460 


0.8920 


0.4687 


0.6490 


0.4948 


0.5342 


3 


0.0692 


1.9706 


0.1446 


1.2283 


0.2219 


0.9285 


4 






0.0307 


2.3249 


0.0674 


1.1638 


5 






0.0094 


4.4003 


0.0188 


2.8050 


6 










0.0039 


4.8755 


7 










0.0018 


8.4742 



3. The attractive Coulomb-interaction integral 



|r-Rfe| 



(b) 



(12) 



where 



— 2Sab{Rij) 



erf(x) 



erf((r^ + r^)i?^,^.) 



■Dab 



X 



e"* dt 



(13) 



is the error function and 



Dab — 



1 + {TJT^)\ 

Using the above formulas, the ground state energy of hydrogen atom takes the form 

{^i\Hi\^i) = 



(14) 



(15) 



a, 6 



^ 3-2V^(rj,)^/^ 2V^TaT,f^ ' 



Minimizing (15) under the condition of proper normalization, i.e. 



„ h y- a b/ 



(16) 



we obtain the values of the parameters Pa and Fq with a = 1, 7. Those parameters for the Slater Is function are 
Usted in Table II. 

The ground state energy of H atom in these three representations arc respectively -0.991 Ry, -0.99912 Ry and 
-0.99987 Ry. Additionally, the expressions (11) and (12), can be used to evaluate the hopping integral < i\Hi\j >. 

The main reason for introducing the Gaussian is to be able to calculate efficiently the many-site integrals involving 
the repulsive Coulomb interactions, which appear in evaluation of U and K parameters of the parameterized Hamil- 
tonian (1) in Sec. 2. In general the Coulomb integral < ij\V\kl >— Vijki has the following form in the Gaussian 
representation 



abed ^ / Ja)Jb) 



yi]ir={9i '9 



V 



(c) (d) 

9 k 9i 



(17) 
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Sac{Rik)SbdiRji)-^ erf {blhcdR-ijki /Babcd) 



where 

Tjabcd — I TD 




i + (rjr,)V V ' i + in/Tdf 



-^2 I -r2\ 1/2 /-r2 I -r2 \ 1/2 



(18) 



&a.c.= (^^) r-^] , (19) 



/p2 I p2 I T.2 I r2\ 1/2 

Babcd = [1±±1a±}^±1a^ . (20) 

Using those expression, the foUowing expressions for parameters U' and K' (i.e. those calculated for the atomic wave 
functions; U and K arc those calculated in the Wannier representation^: 



a,b.c,d 

and 



f/' = = (21) 

„t:^/ '^^/^ V(r^ + rD(rg + r2); b,, 



ihcd 



(22) 



« ^ PaPbPcf3d^eri{bl,,^R,,/ Babcd). 

a,b,c,d 



Explicit calculation of the above integrals requires a number of numerical operations. For example evaluation of 
integral (22) involves computation of 81 integrals in ST0-3G basis and 2401 integrals in ST0-7G basis. Nevertheless, 
they can be calculated efficiently. Obviously, the whole procedure still requires the minimization of the total state 
energy with respect to a in the correlated state. 
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